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PACS 12 . 20 . Fv - Quantum electrodynamics: Experimental tests 
PACS 14 . 80 . - j - Other particles (including hypothetical) 

Abstract. - We investigate the use of ground-based gravitational-wave interferometers for stud- 
ies of the strong- field domain of QED. Interferometric measurements of phase velocity shifts in- 
duced by quantum fluctuations in magnetic fields can become a sensitive probe for nonlinear 
self-interactions among macroscopic electromagnetic fields. We identify pulsed magnets as a suit- 
able strong-field source, since their pulse frequency can be matched perfectly with the domain 
of highest sensitivity of gravitational-wave interferometers. If these interferometers reach their 
future sensitivity goals, not only strong-field QED phenomena can be discovered but also further 
parameter space of hypothetical hidden-sector particles will be accessible. 



Oh Introduction. — Charged quantum fluctuations as 
^predicted by quantum electrodynamics (QED) induce 

onlincar self- interactions of electromagnetic fields [1]. 
This fundamental violation of the superposition princi- 
ple of classical electrodynamics has not yet been observed 
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\^^on the level of macroscopic electromagnetic fields. Even 
though light-by-light interactions have been verified in ex- 

ON|periments involving high-energy photons [2], an investiga- 

^^tion of nonlinear interactions of macroscopic fields would 
probe QED and its vacuum structure in a large-amplitude 
regime which is comparatively little explored in quantum 
field theory. In fact, large-amplitude or strong-field exper- 

T\ iments not only give access to unprecedented fundamental 
^ tests of QED, but also facilitate a search for hypothetical 
particles with light masses and weak couplings to photons 
(hidden-sector searches). This prospect has recently trig- 
gered a remarkable growth of experimental and theoretical 
activities concerned with strong-field and optical set-ups, 
for recent reviews see [3,4]. 

A sensitive probe for vacuum nonlinearities is light prop- 
agation in strong electromagnetic fields. The lowest-order 
nonlinear modifications of Maxwell's theory as induced by 
QED vacuum polarization are described by the (lowest- 
order) Heisenberg-Euler Lagrangian [1], 
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where the scale of nonlinearities is set by the electron mass 
m, a ~ 1/137 denotes the fine-structure constant, and we 
use h = c = 1. The resulting field equations predict that 



a plane wave in a magnetic field B in vacuum propagates 
at a reduced phase and group velocity [5,6], 
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i = -L, II, ay = 14, a_L = 8, 

(2) 

where 9 denotes the angle between the B field and the 
propagation direction. There are two propagation eigen- 
modes polarized parallel || or perpendicular _L to the plane 
spanned by B and the propagation direction. As the 
numerical coeSicients a||^j_ differ for the two polarization 
modes, the magnetized quantum vacuum is birefringent . 
A number of experiments have already been carried out 
[7-11] or designed [12,13] to look for vacuum birefringence 
in terms of high-sensitivity polarimetry. The sensitivity 
limits achieved so far are roughly four orders of magni- 
tude above those necessary for the QED effect. 

An alternative to polarimetry is given by absolute 
phase velocity measurements by interferometry using, e.g., 
gravitational-wave interferometers, as first suggested by 
[14], see also [15]. As gravitational- wave interferometers 
reach their highest sensitivity for optical-path variations 
at a frequency / ^ lO^Hz, the requirements for large-scale 
magnet systems are enormous, as has recently been dis- 
cussed in detail in a concrete proposal in [16]. 

Indeed, using gravitational-wave interferometers for the 
search for nonlinear vacuum phenomena involves many 
parameters which need to be taken into account in an 
optimized fashion. In this work, we propose the use of 
pulsed magnets, as they are developed and used in a sta- 
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ble manner in a number of laboratories worldwide at field 
strengths approaching 100 Tesla. We derive the strain of 
the optical path induced by repeated magnet pulses, as 
this is the central observable at gravitational-wave inter- 
ferometers. We also compute the signal-to-noise ratio for 
selected interferometers as a criterion for measurability. 
As we will demonstrate in the following, not only strong- 
field QED phenomena can be discovered in such set-ups 
but also further parts of the parameter space of hypothet- 
ical hidden-sector particles will be accessible. 

QED-induced strain in gravitational-wave inter- 
ferometers. — Gravitational wave interferometers con- 
sist of two perpendicular interferometer arms of equal 
length L. Mirrors at the arms' ends form an evacuated 
cavity for two in-phase laser beams one for each cavity. 
An incoming perturbation such as a gravitational wave 
leads to a relative change AL of the length of the interfer- 
ometer arms, which manifests itself as a phase difference 
of the two laser beams. The strain h{t) = AL/L then cor- 
responds to the amplitude of the incoming perturbation. 
Remarkably, present day interferometers are, in certain 
frequency ranges, sensitive to a strain of /i ^ 10^^^: fu- 
ture upgrades aim at further improvements of one or two 
orders of magnitude. 

In the context of quantum-induced refractive properties 
of the vacuum, it is not the arm length but the corre- 
sponding optical path length L which can be modified by 
a strong field. If an external field is applied in a region 
X < L in one interferometer arm, an optical-path differ- 
ence AL is induced, (AL/x) = 1 — v. Prom now on, we 
concentrate on the parallel mode for which the velocity 
shift is maximal. For this mode, the resulting strain is 

. |(.3x:o-^.,(M)\ 

For a sizable strain, the magnet-length-to-detector-arm ra- 
tio x/L should be maximized. As x is constrained by the 
technical realizability of the magnetic field coil and the 
detector sensitivity requires long arm lengths, a suitable 
compromise between detector arm length and sensitivity 
has to be found. 

Detector sensitivity. — The detector sensitivity to 
the relative shift of the length of the interferometer arms 

AL/L is limited by various sources of noise. At low fre- 
quencies / < 40Hz, the main limitation arises from seis- 
mic activities for ground-based detectors. For instance 
at / ^ lOHz. which seems ambitious but feasible for big 
dipole magnets, the sensitivity measure Sh{f) (as detailed 
below) is suppressed by more than three orders of magni- 
tude compared to the peak sensitivity and depletes rapidly 
for even smaller frequencies. At intermediate and higher 
frequencies, thermal and shot noise, respectively, limit the 
detector sensitivity [17,18]. We conclude that observing 



the quantum-induced phenomena requires a magnetic field 
modulation in the optimal frequency range for a given de- 
tector. Pulsed magnets which c;an be modulated over an 
appropriate parameter range therefore appear to be good 
candidates. 

Information about the optimal frequency range is en- 
coded in the total spectral density function Sh{f) of the 
noise, see e.g. [18]. For our estimates, we concentrate 
on advanced LIGO (in preparation) and GEO600 (op- 
erational). The projected sensitivity of both detectors 
slightly depends on the details of the event acting as a 
source for the interferometric signal. In the case of LIGO, 
we use the strain sensitivity for neutron-star binaries [19] 
for our calculations, which is satisfactory over a wide range 
of intermediate frequencies (for a typical expected sensi- 
tivity curve of the advanced LIGO, see [21]). For GEO600, 
we use the typical sensitivity data available at [20] , which 
can be well approximated by a fit function of the form 

Sh{f) = Soiifo/fr + 2(///o)^= + 2)/5, (4) 

where fo = 560Hz, 5*0 = 7 x IQ-'^Rz-^, pi = 3.8, P2 = 3 
near the sensitivity maximum. For the advanced LIGO, 
no simple fit is available. As the signal induced by QED 
or other hypothetical particles can be predicted exactly, 
we expect that the interferometer sensitivity can even be 
optimized for the present phenomenon. 

As a measure for the observability of a shift of the op- 
tical path AL/L, we determine the signal-to-noise ratio 
(SNR) d of the induced strain. In the present case, the 
SNR equals the expectation value of the detector output 
divided by the standard deviation of the output variable 
due to noise. Using a matched filter (or "Wiener filter") 
for the signal, the SNR d for a gravitational wave inter- 
ferometer is given by (see e.g. [18, 22, 23] and references 
therein) 

where h{f) is the Fourier transform of the induced strain. 

Pulsed magnetic fields. — As an example, we con- 
sider pulsed fields that can be obtained at the Dresden 
High Magnetic Field Laboratory (HLD) [24]. The HLD 
aims at providing lOOT fields generated by a solenoid in 
a non-destructive set-up, i.e., the infrastructure is main- 
tained and the experiment can in principle be repeated 
arbitrarily often. 

As the magnetic pressure is given by Pma.g = -B^/2/xo) 
already at fields strengths of about B = SOT, the pres- 
sure on the coils is four orders of magnitude above the 
atmospheric pressure, which requires a careful coil design. 
As a consequence, the coils are usually heavily mantled 
and it is difficult to render the interferometer laser beam 
orthogonal to the external magnetic field within a single 
coil setup. In order to maximize the shift of the optical 
path, a pair of Helmholtz coils has to be used instead. 
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In this manner the laser beam of the interferometer can 
be aligned in parallel to the field coils and thus mainly 
orthogonal to the magnetic field lines without interfering 
with the coil mantle. 

For such "spUt coils", a maximum field strength of up 
to 60T appears technically feasible, at a coil diameter of 
about X = 0.2m and a coil separation of O(lcm). As 
the beam waist of the interferometer lasers is of the or- 
der of cm, the interferometer beam can fit between the 
magnet coils, even though the issue of stray photons may 
require further discussion which goes beyond the scope of 
the present work. Also for standard Helmholtz setups, the 
field is roughly constant at a sizable extent only along the 
direction of the magnetic field lines, whereas the quantum- 
vacuum effect requires a sizable field length perpendicular 
to the field lines; the length of the latter is of the or- 
der of the coil separation. For the proposed setup, the 
coil design thus needs to be optimized to provide for high 
(but not necessarily constant) magnetic field strengths, 
spatially extending orthogonally to the direction of the 
field lines. 

A typical pulse undergoes a damped oscillation with 
pulse frequency u-q and damping rate 7. For A'' subse- 
quent pulses at times Iq. . . t^-i, a satisfactory description 
is given by 



B{t) = 5, 
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ti) sin(27ri^B(< — ti)) exp{—j{t — ti)). 

(6) 

Here, we have ignored that successive pulses have no tem- 
poral overlap in a single-magnet set-up. However, as the 
pulse repetition rate vp = — ti) is much smaller 

than the damping rate (see below), Eq. (6) is a well 
justified approximation. The pulse frequency i/b in Eq. 
(6) depends on the total capacity of the capacitor banks 
and can lie in the range ©(ms . . . s), while the damping 
rate 7 is mainly determined by the heat capacity of the 
coil. In addition, the achievable pulse repetition rate fp 
in a non-destructive mode depends strongly on the desired 
peak field strength. An ambitious, but nevertheless feasi- 
ble Helmholtz setup should be able to achieve a maximum 
field strength of Bmax = 60T, followed by a reverse field 
of i?min = — 6T and thus a damping to about 10% of the 
peak field strength ^. This choice fixes the amplitude of 
the model pulse (Eq. (6)) Bq 148T and implies the 
constraint 

o 1 

7 = zvb In 



(7) 



We use the remaining free parameter vb for optimizing 
the SNR (5) (within the technical limitations). For this, 
we need the modulus of the Fourier transform of the 
strain. Ignoring overlap terms of successive pulses as ar- 



^It can be expected that also pulses with damping to about 70% 
will be obtainable in the near future. The above pulse parameters 
are in that sense conservative, since a smaller damping factor 7 leads 
to a higher strain (cf. (8) 
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gued above, we find 
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For a single pulse N = 1, the trigonometric functions in 
Eq. (8) cancel and the dependence on the repetition rate 

i/p drops out, as expected. For a large number of pulses 
N, the trigonometric functions yield a representation of a 
6 comb. 
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At large N, only frequencies which are multiples of the 
pulse repetition rate //b thus contribute to the SNR (5). 
As ub is much smaller than the frequencies dominating the 
SNR, the contributing frequencies form a quasi-continuum 
such that the sum in Eq. (9) c;an well be approximated 
by an integral. As a result, the square of the SNR for A'' 
pulses can be expressed in terms of the single-pulse result 
to a good accuracy: 
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The reproducibility of the signal by non-destructive pulsed 
magnets thus is a lever arm for an enhancement of the SNR 
by a factor of V^. For the advanced LIGO with L = 
4000m, the sensitivity curve has a broad minimum of the 
order of Shif) ~ 10~^^Hz~^ for frequencies ranging from 
approximately 50Hz to 500Hz. Maximizing (P by varying 
the pulse parameter ub yields ub ~ 47Hz, implying 7 « 
217Hz by means of Eq. (7). Inserting these values into Eq. 
(5), we obtain the SNR for a single pulse. 
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As a result, about N w 2763 pulses are required in order 
to achieve a total SNR of 0(1). Depending on the de- 
tails of the setup, an SNR of 0(10) might eventually be 
required. For the following feasibility study, however, we 
only demand for an SNR of 0(1). This is also justified 
because the expected signal can be predicted to a high 
accuracy which will allow for an adapted noise filtering. 

As mentioned above, the re-cooling time for the magnet 
system which determines the pulse-repetition rate depends 
mainly on the pulse energy. A realistic estimate lies in the 
order of several minutes, implying a continuous operation 
of the facility for a few days. This appears perfectly feasi- 
ble. GEO 600 is considerably less sensitive than advanced 
LIGO but, for our purposes, profits from the shorter arm 
length of 600m. Maximizing (P with respect to ub yields 
Ub ~ 273Hz with 7 « 1259Hz and thus a pulse length be- 
low 1ms. As a result, « 2.3 x 10^ pulses are necessary 
to observe the QED induced strain at GEO. This corre- 
sponds to an unrealistic measurement time of a few years. 
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In consequence, GEO in combination with presently avail- 
able pulsed magnets is not well suited for the observation 
of the QED induced strain. Nevertheless, it still has a 
new-physics discovery potential, see below. 

Search for hidden-sector particles. The investi- 
gation of vacuum nonlinearities can also be used to search 
for hypothetical particles with small masses and weak cou- 
plings to photons [4]. An apparatus designed to observe 
the QED effect will also explore a significant range of a 
new-physics parameter space. Here, we concentrate on 
potential velocity shifts induced by axion-like particles 
(ALPs) [25] or minicharged particles (MCPs) [26]. 

MCPs with charge Q = ee and mass rrie induce vac- 
uum eS'ects analogous to electrons. As the MCP mass 
can be much lighter than the electron, quantum-induced 
velocity shifts have to be calculated to all orders in the 
field strength parameter eeB/rn^ and the frequency w/m^ 
[27, 28] . MCP masses corresponding to a Compton wave- 
length larger than the volume inside the Helmholtz coil 
cannot be fully resolved within the setup. Thus, the 
coil separation of ©(Icm) constrains the search for MCP 
masses to ^ 2 x lO^'^eV. Using the formulas of [26,28] 
for a Dirac spinor MCP (scalar MCPs behave similarly 
[28]), we calculate the effect for the polarization compo- 
nent parallel to the external field, which maximizes the 
velocity shift as in the QED case. For an interferometer 
laser with uj = 1.2eV and the pulse shape as used for the 
QED effect, we obtain exclusion limits in the fractional 
charge-mass plane {e^rrie) by demanding a SNR of 0(1), 
see Fig.l. These exclusion limits display two characteris- 
tic limits which correspond to the two asymptotic limits 
of the velocity shift: \ — v^ e^B^/m'^ for large masses (cf. 
Eq. (2)) and 1 - w e^/^B^/^/w''/^ for small masses. 

We find that already a single-pulse {N ~ 1) measure- 
ment at advanced LIGO can approach the current best 
laboratory bounds [28] derived from PVLAS data [8] with 
potential minor improvements in the larger-mass range. 
Assuming a measurement time of ten days at GEO with a 
magnet re-cooling time of 5 minutes, implying N ~ 2880 
pulses would arrive at a similar bound. The small-mass 
asymptotics of the GEO bounds for iV = 1 and N = 2880 
are also shown in Fig. 1. Advanced LIGO with N = 2763 
as needed for QED can reach a sensitivity of e ~ 10^^ for 
nie < O.OleV. This would compete with current cosmolog- 
ical bounds [29]. Of course, astrophysical energy loss con- 
siderations can lead to much stronger bounds [30] reaching 
down to £ ~ 10~^*. However, as stellar physics probes a 
much larger momentum-transfer regime, these bounds de- 
pend on the underlying microscopic physics and thus may 
not apply to laboratory experiments [31]. We conclude 
that an investigation of the QED effect can also improve 
current laboratory bounds on MCPs to an unprecedented 
level. 

Axion-like particles (ALPs) are uncharged scalars (S) or 
pseudo-scalars (P) with a mass and an effective ALP- 
photon interaction of the form Cp = g(j)pE ■ B or Cg = 
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Fig. 1: New-physics discovery potential for spin-| minicharged 
particles (MCP) in the coupling-mass plane. For MCPs, al- 
ready a single-pulse (A'' — 1) measurement at advanced LIGO 
(or A*' ~ 2880 at GEO) can approach or slightly improve 
the current best laboratory bounds from PVLAS. The use of 
A'^ = 2763 pulses as needed for the QED effect can lead to 
sizable improvements over the whole mass range. 



|f;0s(B^ — E^) with coupling g. The interaction structure 
implies that scalars only couple to the _L mode, whereas 
pseudo-scalars couple to the || mode of the propagating 
light in a magnetic field. Loosely speaking, an effective 
reduction of the phase velocity arises from the fact that the 
corresponding photon partly propagates as a massive ALP 
component. The corresponding velocity shifts read [25] 



1 - = 1 - «i = 
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where oj denotes the laser frequency, and x is the length 
of the magnet interaction region. As discussed in detail 
in [28], the polarimetry of strong-field particle searches 
can distinguish between the various particle scenarios, see 
also [16]. For the current study, we only concentrate on the 
velocity shift irrespective of the polarization dependence 
for simplicity. 

Let us first consider the ALP parameter range for the 
coupling g and mass that can be probed by a single 
pulse. Using the pulse shape as for the QED case, we 
obtain the accessible region in the mass-coupling plane, 
see Fig. 2 for a comparison with the current best labora- 
tory limits from PVLAS and GammeV [32] (see also the 
BMV experiment [33]). Already a single-pulse measure- 
ment at advanced LIGO can improve on existing bounds 
for masses above ImeV. For ^ 4 x 10~^eV, higher 
number of pulses as needed for the QED effect give access 
to a parameter space which is significantly larger than cur- 
rent laboratory limits. Again, astrophysical considerations 
and observations lead to stronger bounds by a few orders 
of magnitude in the g ~ lO^^^GeV"^ regime [34], but may 
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Fig. 2: This figure shows the discovery potential in the 
coupling-mass plane for axion-like particles (ALP) at LIGO 
and GEO, respectively. For ALPs, the current best laboratory 
bounds by GammeV and PVLAS can be improved significantly 
in the meV mass range and above by advanced LIGO as well 
as already by GEO. 



der of magnitude in the field strength and lose an order of 
magnitude in the field extent, the quantity xB^ can gener- 
ically still be an order of magnitude larger for pulsed fields 
than for dipoles. A similar consideration has inspired the 
development and use of pulsed magnets in the BMV ex- 
periment [10, 33] which finally aims at a parameter goal 
of xB^ ~ 600T^m (recent experimental results of BMV 
have been achieved with xB^ ~ 40T^m) [10]. The pulsed 
Helmholtz coil configuration considered in this work which 
is inspired by ongoing experiments at the Dresden High 
Magnetic Field Laboratory would yield xB^ ~ 720T^m, 
which agrees with the design goal of BMV. But whereas 
the time dependence of the field pulse at BMV can be a 
disturbing factor for the BMV polarimetry, the time struc- 
ture of the field pulse is a necessary and advantageous key 
feature for a strong-field experiment at a gravitational- 
wave interferometer. 

Given the prospect of exploring a new parameter regime 
of strong-field quantum field theory with implications for 
the search for new elementary particles, establishing a 
strong-field quantum-vacuum program at gravitational- 
wave interferometers appears to be worthwhile. 



not be directly applicable to laboratory experiments due 
to a different momentum-transfer regime [31]. 

Conclusion. — The sensitivity goal of ground-based 
gravitational-wave interferometers appears well suited to 
use these experiments for exploring the strong-field do- 
main of QED. In order to generate a strongly magne- 
tized quantum vacuum in such an interferometric exper- 
iment, we have identified pulsed magnets as an advanta- 
geous strong-field source for two reasons: they provide ex- 
tremely strong laboratory magnetic fields, and their pulse 
frequency can be perfectly matched with the region of 
highest sensitivity of the gravitational-wave interferome- 
ters. 

For our quantitative estimates, we have concentrated 
on the advanced LIGO detector, as its sensitivity goal 
matches with currently available field strengths already 
in a rather conservative estimate. Pushing the various 
components to their limits may facilitate a detection also 
at the gravitational-wave interferometers, which are cur- 
rently operational such as GEO 600. Also the fact, that 
the quantum-induced signal can be theoretically predicted 
to a good accuracy may give rise to an improved noise fil- 
tering. 

From a general perspective, the QED velocity shift as 
well as the MCP signal in the large-mass domain and the 
dispersive ALP effect scale with xB^, where B is the am- 
plitude and X the extent of the magnetic field. For the use 
of gravitational-wave interferometers, also a suitable time 
variation of the magnetic field is needed. Whereas pulsed 
fields profit from extremely high fields and a suitable time 
variation, their deficit is a smaller extent in comparison 
to dipole magnets. Since pulsed fields win roughly an or- 
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